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Abstract 

We reassess a work done by Migdal on the spectrum of low-energy vector mesons in 
QCD in the light of the AdS-QCD correspondence. Recently, a tantalizing parallelism was 
suggested between Migdal's work and a family of holographic duals of QCD. Despite the 
intriguing similarities, both approaches face a major drawback: the spectrum is in conflict 
with well-tested Regge scaling. However, it has recently been shown that holographic 
duals can be modified to accomodate Regge behavior. Therefore, it is interesting to 
understand whether Regge behavior can also be achieved in Migdal's approach. In this 
paper we investigate this issue. We find that Migdal's approach, which is based on 
a modified Pade approximant, is closely related to the issue of quark- hadron duality 
breakdown in QCD. 
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1 Introduction 



The recent years have seen an emergence of 5-dimensional modelling in an effort to make con- 
tact with QCD, the main motivation under most of the models relying upon the gauge/gravity 
duality conjectured by Maldacena almost a decade agoU 

Even though the original Maldacena conjecture related a slice of superstring theory on 
anti-de Sitter (AdS) background with M = 4 SYM, both the belief that QCD is rooted in 
string theory and, especially, the high-energy conformal symmetry of QCD yield room for 
hoping that the conjecture may be extended beyond its actual formulation to provide new 
insights into the strong interactions. Indeed, rather simple models capturing some of the 
features of QCD like chiral symmetry breaking and confinement are in remarkable agreement 
with phenomenology. 

The common setting of these models is a 5-dimensional AdS background. A 4-dimensional 
boundary brane is introduced, upon which Standard Model fields are bound to propagate. 
Confinement is modelled by cutting the bulk space down at another boundary brane, thus 
introducing an infrared scale that mimics Aqcd- Spontaneous chiral symmetry breaking 
[3j [U [5], quark masses [U [6] and OPE condensates [7] have also been implemented over 
the last yearsEI However, the manifestly wrong scaling of the predicted spectrum, namely 
m 2 ~ n 2 , has been a burden to all these models, not only from the phenomenological side 
but also because their eventual connection to string theory was doubtful. 

Only recently was it realized [10] that the wrong mass-scaling was not an essential draw- 
back of holographic models of QCD, but due to the choice of boundary conditions. One can 
draw an analogy between the 5-dimensional Kaluza-Klein modes in the bulk of the afore- 
mentioned models and the quantum mechanical excitations in a square well potential. By 
changing the profile of the boundaries one can dial the energy spectrum at will and should be 
able to get to the desired mass-scaling. In the language of holographic QCD, this amounts 
to a dynamical dilaton, scaling (asymptotically) as &(z) ~ z 2 , whereas the metric is kept 
(asymptotically) purely AdS. Therefore, hard- wall models, i.e., those with a sharp infrared 
brane, yield m 2 ~ ra 2 , as opposed to infrared-improved models, where the asymptotic form 
of the dilaton ensures m 2 ~ n. 

It turns out, as first pointed out in [11], that hard- wall holographic models bear a tan- 
talizing resemblance with a work done by Migdal [12J almost 30 years ago on modelling the 
spectrum of vector mesons in QCD0 This very ambitious project relied on linking short 
distance QCD to the meson spectrum through quark-hadron duality. The information on 
short distances was used to make a Pade approximant to the whole correlator and thereby 
information on the spectrum was extracted. This procedure was designed to be a systematic 
step by step approximation, in the sense that adding more information on short distances 
one should get closer to the true QCD spectrum. Thus, at least formally, the programme 
was improvable. Therefore, when the Pade for the leading order perturbation theory yielded 
a spectrum with poles at zeros of Bessel functions, meaning that the tower of vector mesons 
in QCD should scale as m 2 ~ n 2 , the apparent contradiction with experimental evidence was 
seen as an artifact of the truncation. Should more orders in the OPE be added, the spectrum 
would smoothly reshuffle to display Regge trajectories. 

The similarity between Migdal's result and that of hard-wall holographic models is, in 

2 See, for instance, Q3H] f° r comprehensive reviews. 

3 See also, e.g., [HE] for the effects of modifying the AdS metric in the infrared. 

4 Recently it was also claimed [13] that the extension of Migdal's programme to three point functions carried 
out in Q3] has its counterpart in the holographic light-front wavefunctions of |15| . 



1 



a sense, not so surprising. In the simplest hard-wall models the 5-dimensional background 
is filled with a pure AdS metric, meaning that the correlator looks like that of a free quark 
theory. Both approaches therefore start from the same conformally invariant correlator. The 
interesting thing is that the analogy seems to go further and one might even argue |16j that 
holographic models can be interpreted as the limit of a series of Pade approximants. A 
qualitative understanding of this should be found in the deconstructed version of holographic 
QCD models [17] , where an open moose model of hidden local symmetries [18] was shown to 
be interpretable as the latticized version of a fifth dimension. A [N, N] Pade approximant 
would then be the dual picture to an open moose with N hidden symmetries. 

In [16], Migdal's procedure was intended to shed some light on AdS-QCD. Here we reverse 
the point: we reassess Migdal's work with the benefit of using the geometrical intuition 
of AdS-QCD models. A natural question then arises: if Regge scaling was achieved in 
holographic models through a non-trivial dilaton background, what is the analog in Migdal's 
approach? We will see that one is forced to include non-perturbative effects beyond the OPE, 
generically referred to as quark-hadron duality violations in the QCD literature [19J. Actually, 
we will show that non-perturbative effects were already present in Migdal's original work. 

The discussion on quark-hadron duality breakdown will also allow us to grasp the meaning 
of the dilaton's quadratic scaling in holography, an interesting question partially answered 
in the context of string theory [TO]. Still, one would like to have some physical intuition 
in the context of QCD. We explicitly show that the dilaton's infrared profile in Regge-like 
holographic models is dual to a specific modeling of duality breakdown in QCD. 

The outline of the paper is as follows: sections 2, 3 and 4 deal with 4-dimensional models: 
in section 2 we introduce a toy model with explicit Regge scaling, while in section 3 we discuss 
Pade approximants. Section 4 is devoted to an analysis of Migdal's modification of the Pade 
approximant and the relationship with the issue of quark-hadron duality violations. We will 
work in the large- A^ c limit of QCD. There, resonances are narrow and therefore mass-poles are 
well-defined. The aim is to use our toy model, where results can be worked out analytically, 
as a probe to reassess the accuracy of Migdal's programme. Contact with 5-dimensional 
holographic models is made in the last section. 

2 A toy model 

Following [12], we consider the following two-point vector correlator 

U^(q) = i J d 4 x e ^(0|T{V^(x)(V^)t (0)} | 0) f (1) 

with currents defined as 

V^(x) = u{x)rd(x) . (2) 
In the chiral limit, Lorentz symmetry implies that 

K l '(l) = (Q tl Q V -1 2 9n^v(q 2 ), (3) 

namely that there is only one form factor, which satisfies the following dispersion relation 

f°° dt 1 

(q 2 ) = q 2 / — — -r ~Im U v (t) + U v (0) , (4) 

Jo t(t - q l - le) vr 
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where Ily(O) is a subtraction constant. In the strict large- N c limit, a two-point Green's 
function consists of single particle exchange of an infinite number of stable mesons [20] . 
Therefore, the spectral function reads 

1 °° 

-Im U v (t) = Y, F n*(t - M 2 (n)) . (5) 

n=0 

Unfortunately, this is the furthest one can go in full generality. Decay constants and masses 
in large--/V c QCD are unknown. However, phenomenology in N c = 3 QCD seems to indicate 
that Regge behavior holds and, in the vectorial channel, already sets in from the very first 
resonance. Therefore, we will choose as an educated ansatz |19j 

F 2 = F 2 , M^(n) = m 2 v + an . (6) 

Inserting ([6]) back into we end up with the following expression for the correlator 



2 



iMs 2 ) = - 

a 



(7) 



where is the digamma function, defined as^| 



^(0 = ^iogr(0- (9) 

The free parameters F and a can now be fixed upon matching with the first terms in the 

OPE for Uv(Q 2 = —Q 2 )- Indeed, in the Euclidean half-plane, the digamma function admits 
the following integral representation 

r°° ( e~ l \ 

* (() = i dt {--T^)' (10) 

from which the large-£ behavior can be extracted 



00 R 

£>■ 



2/i 



2£ ^2nf 2n 

n=l 

Thus, the OPE for our model reads 



nTO = 4,„ g ^ + |-, (12) 



where the condensates are given by [21 



C2k = (-l) k ^-B k (^) . (13) 



Notice the connection between this model and the Veneziano model (Beta function) for fixed spin, i.e., 

B(n ' s) - T(n + s) ■ (8) 
In both cases, the mass-poles are dictated by the poles of the Gamma function. 
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and -Bfc(£) stand for the Bernoulli polynomials. On the other hand, in real QCD, to leading 
order in a s , 

n^(Q 2 ) = -^iog^ + ... (14) 



Thus, we immediately get 

4 N c F 2 
3 (4tt) z a 

As it stands, the model only fullfils one short-distance constraint. In principle, one could 
refine the model in order to match the first OPE condensates. This can be easily achieved 
by relaxing our initial ansatz for the decay constants ([6]) and letting an arbitrary number 
of them as free parameters. For such extensions of our model, we refer the reader to |21| . 
However, in this paper we will not be interested in the phenomenological implications of the 
toy model. For all our considerations we can safely stick to the simplest version of Eq.([7]). 

3 Building the Pade approximant 

With our model for the two-point vector current correlator of ([7]), we can now attempt to 
build its Pade approximantjj With the advantage of knowing the final answer, we will try 



to assess the validity of the procedure followed in [12]. We begin with some notation and 
definitions. 

A [N, M] Pade approximant to a function is defined as the quotient of two polynomials 
Vm and Q N 

such that the first iV + M + 1 derivatives of the function at a point match those of the 
approximant. With full generality, we can set Qtv(O) = 1. Notice that one should be able 
to define (at least formally) a Taylor expansion of the function around a point, even if the 
expansion is divergent. The next issue to address is that of convergence. One would wish 
that the [N, M] rational should converge to the function as N,M — > oo. This is true for a 
class of functions known as Stieltjes functions. Stieltjes functions admit the following power 
expansion about the origin 

oo 

f(z) = J2fnZ n , (17) 
n=0 

with the coefficients (also called moments) given by the following integrals 

j_ 

/n = (-l)" l R U n v(u)du, (18) 





where R is the radius of convergence of (|17p and v{u) is some positive definite integration 
measure. Therefore, ()1T|) admits the following representation as a Stieltjes integral 

/(z) = /o TT^" (t,) - (19> 



6 Pade approximants have been a useful tool for widespread applications in condensed matter and particle 
physics. See, for instance, the very interesting examples of [22] . 
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With the change of variable u = t 1 one obtains 

f°° dt 
l R t(t + z) 



which matches (jlj) with the following identifications^ 



' 1 1 



-imn v (t). 

7T 



(20) 



(21) 



Therefore, as recently pointed out in [23], the vector vacuum polarization admits a represen- 
tation as a Stieltjes integral. This is most fortunate, and a series of interesting consequences 
ensued First and foremost, it ensures convergence of the Pade approximant of (]16p for the 
special subset \{M + J),M], J > —1. The convergence region covers the whole complex q 2 
plane, save for the Minkowski axis, where the poles sit. Furthermore, all the poles of the 
approximant are simple and located on the negative real axis, i.e., on the physical axis, and 
residues are positive. 

Therefore, from a physical point of view, the Pade approximant to ([7]) can be built, at 
least formally, out of the infinite number of terms of the chiral expansion. Pade poles can then 
be interpreted as meson mass-poles whereas Pade residues play the role of the meson decay 
constants. Convergence then means that our Pade approximant has to match asymptotically 
the real QCD spectrum. 

The strategy followed in |12| was to construct the symmetric [N, N] approximant around 



an expansion point q = —[i on the far Euclidean axis. The Pade equations then read |12lH6 



d(q 



2\n 



Uv(q 2 )QN(q 2 )-V N (q 2 



n 



0,...,2N. 



(22) 



It is worth noting at this point that the convergence theorems above no longer apply if one 
builds the Pade away from the origin. Interesting exceptions are meromorphic functions. 
There it can be proved that convergence is guaranteed as long as the Pade is constructed on 
any non-singular point of the real q 2 axislfl In the strict large- N c limit Green's functions are 
known to be meromorphic, which means, in particular, that with the model of section 2 the 
Pade constructed around an arbitrary q 2 = — [i 2 converges. However, for finite N c nothing 
compels our Pade to approach the original function. It is true that often Pade approximants 
converge even when there is no theorem to guarantee their convergence, but this pragmatic 
approach is far from being rigorous. 

Here we are just interested in the mass spectrum, which is encoded in Qn(q 2 )- We will 
therefore keep only the upper half of the previous set of equations, where the contribution of 
Vn(q 2 ) vanishes, namely 

"'" " ' ' ^ ' 0, n = N+l,. 



d{q 2 )n n v (q 2 )QN( q 2 ) 

Use of Cauchy's integral formula leads to 
d 



d(q 



2\n 



n V (q 2 )QN(q 2 



2~ni 



dq 2 



(q 2 + v 



2\n+l 



,2N 



U v (q 2 )Q N (q 2 



(23) 



(24) 



7 Notice that our ansatz for the spectrum ([5]) assumes a mass gap with the first resonance sitting at q 2 



Therefore, R = my in (|20[1 . The spectral integral beginning at t = 
We will hereafter use the former notation for simplicity. 
8 See, for instance, [24] . 

9 I thank G. Lopez-Lagomasino for helping me clarify this issue. 



: m v . 

or t = niy thus yields the same result. 
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(? 2 = V) 




Figure 1: Contour integration chosen in {24]) ■ The singularities on the right half-plane are the poles 
of TLv(q 2 ), sitting at q 2 = m 2 , + ak for natural k. The pole on the left half-plane is the pole with 
multiplicity (n + 1) explicitly shown in the denominator of {24] - The circle is defined by \q 2 \ = sq, 
which will be eventually sent to infinity. 



where 7 is the contour depicted in figure 1. We can relate the previous expression to the 
spectral function by using the conventional strategy in sum rules of splitting the contour in 
two paths (see, for instance, |25j). In N c = 3 QCD one can readily show that the integral 
over the circle l^ 2 ) = sq vanishes when the contour is sent to infinity, and only the integral 
over the physical axis contributes. This is related to the fact that the asymptotic behavior 
of the spectral function goes to a constant and sc0 



lim -ImIL/0) Q N (t) ~ t 

t— >oo 7T 



N 



(25) 



Therefore, N + 1 derivatives are sufficient to get a Kallen-Lehmann representation for the 
Pade equations, namely 



d{q 2 ) n 



U v (q 2 )Q N (q 2 



dt 



(* + A* 



2\n+l 



-ImU v (t)Q N (t) 

IT 



(26) 



n 



N+1,...,2N, 



where we made use of Schwarz reflection formula and the reality of Qn(q 2 )- Following 
Migdal's strategy, at that point we will rely on (local) quark-hadron duality and assume that 
the spectral function can be inferred from the OPE of f)12[) . Then, 



-ImU v (t) ~ ^7^2 

7T 3 (4:1T) Z 



: (2fc+2) 



(27) 



k=0 



where 5^ k \t) stands for the A;-th derivative of the Dirac delta and S^°'(t) = S(t). If we just 
keep the constant term in (|27|) -i.e., the logarithmic piece in (|12|) ~. then the Pade equations 
(f26|) reduce to 

-^+rSiv(i)=0, n = JV + l J ...,2JV. (28) 







10 In the strict large-Af c limit, equation (|25|) is ill-defined because of the discrete nature of the spectral 
function. However, one can show that (TZU)) still holds. See the Appendix for a detailed proof. 
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By reorganizing the powers of (t + fj?) in (|28p . Migdal realized that Qn(q 2 ) are proportional 
to Jacobi polynomials. Actually, it seems that the first derivation of the Pade approximant 
to the logarithm was performed by Jacobi, Rouche and Gauss q Their result for the poles 



was 



Q N (q 2 



2F1 



-N,-N:U- 



+ pirn ( a 



N 



/i 2 + q 2 



(29) 



where 2^1(0, b; c; d) is Gauss' hypergeometric function and P^'°^ stands for the (0, 0)- Jacobi 
polynomial, i.e., the Legendre polynomial. Up to an overall constant out front, this last 
expression is the one found in the recent analysis of |16j . Plugging (|29p in (|22p we can solve 
for 7-V(g 2 ). Assembling all the pieces, the final result for the correlator reads 



(q 2 + V 



2\N 



N 

£ 

k=0 



H 



N-k 



2\ k 
2 



X 



» 2 -q 2 
/U 2 + q 2 



(30) 



where are the harmonic numbers. The expression above is the Pade approximant to the 
logarithm, i.e., when N —* 00 the logarithm we started from is recovered. 



4 Migdal's limit and quark-hadron duality 



However, Migdal went one step beyond the strict Pade approximation. At low energies 
(q 2 <C /J?) the argument of the Jacobi polynomial above can be Taylor-expanded to yield 



lim 



[i 2 - q 2 



g 2^ /U 2 + q 2 



2I 

2 



Then, using the asymptotic formula [27] 



lim P- 



(0,6) 



AT— >oo 



N 



2N 2 



MO 



one ends up with 



Q N (q 2 



Jo 27V, 



(31) 



(32) 



(33) 



and the full correlator takes the form 

^ / 2n 4 N r 



3(4vr) 



log 



7T 



Y (qA 



H 2 J (qA ) 



A 



2N 



(34) 



It is interesting to examine the last expression both in the Euclidean and Minkowki regions. 
In the far Euclidean it can be shown that the correlator reduces to (Q 2 = — q 2 ) 



n v (Q 2 ) 



4 N r 



3(4vr) 



O 2 

log^ + 0(e 



-2QA> 



(Q 2 > 0). 



(35) 



The last term does not stem from the Pade approximation. Rather, it is a consequence of 
the combined q 2 <C /x 2 and 7V — > 00 limits taken above. Despite its exponential suppression, 
the new term is of the outmost importance. First and foremost, it does not belong to the 



'See 26 for details. 
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OPE, which is a series expansion in inverse powers of Q 2 . The presence of exponentially 
suppressed terms in the Euclidean half-plane is a reminder that the OPE is not an allowed 
expansion in the whole q 2 complex plane. At least on the physical axis the expansion breaks 
down. Generically one writes 

U v (q 2 ) ~ II« Pi V) + A(<Z 2 ) , (q 2 » 0), (36) 

where A(g 2 ) is referred to as quark-hadron duality violation. This extra term can be under- 
stood on the basis that the OPE, being a regular expansion, loses track of the precise details 
of the spectrum, which are then encoded in the A(q 2 ) function. A(q 2 ) thus contains the 
singularities (cuts and poles) of Hv(q 2 )- Even though the origin of such extra term is well 
understood, its characterization is a most challenging task. 

Two basic patterns have been studied so far in the literature, associated with infinite-size 
and finite-size singularities [19]. In the former case, based on a model of resonances, the 
authors of [211 12H] concluded that 

A(q 2 )^H(q 2 )e-^^ . (37) 

where H(q 2 ) and F(<f>, N c ) are in general unknown. The only generic constraint comes from 
the fact that in the strict large-iV c limit the spectrum is an array of Dirac deltas. This 
requires that J-{4>, N c — > oo) has to cancel on the Minkowki axis (4> = 0, 2tt), therefore killing 
the exponential damping of Eq. ([37|) . 

Finite-size singularities were studied in an instanton model |19j . Quite unexpectedly, the 
construction of Migdal, even though it is claimed to be a model of resonances in the large- N c 
limit, leads to an instanton-like pattern for duality violations, as first noted in [7] 1^1 

Obviously, the pattern is induced by the appearance of the infrared regulator A, which 
arises after the limits q 2 <C fJ 2 , N — > oo have been taken. It is important to notice that 
both limits are correlated, i.e., fi/N is to stay constant. This peculiar continuum limit spoils 
the convergence of the Pade approximant. The best way to visualize this is to consider the 
expression for the Pade residue 

nv(q 2 )QN(q 2 )-V N (q 2 )=K N (q 2 ), (38) 

in terms of which our correlator can be expressed as 

u i 2, ?W) 4 N c -q 2 Tl N {q 2 ) 

QjvOr) 3 (4vr)^ n z QN{q 2 ) 

For the logarithm it can be shown that the last term above is a Meijer G function which 
vanishes in the continuum limit (N — > oo), thus ensuring that our Pade approximant is 
indeed convergent. Notice, however, that fixing the ratio [i/N gives a non-vanishing residue, 
cf. Eq. (|34p , Convergence of the Pade approximant is lost due to a source of duality violation 
put by hand and parametrized in terms of an infrared scale A. 

Let us now turn our attention to the Minkowski axis to analize the impact of the new 
scale A. According to our previous discussion about the duality violating pieces, they have 
to account for the singularities of the correlator. This is indeed what we observe for the toy 
model of section 2. The asymptotic expansion of Eq. (|llh leading to the OPE of Eq. (112j) is 



12 It is worth noting that instanton-like duality violations are expected to be suppressed at large- N c , as 
Migdal himself noted. Duality breakdown due to resonances, however, is a leading order effect in the 1/N C 
expansion. 



s 



no longer valid in the physical axis, where the digamma function has its poles. This can be 
circumvented by recalling the reflection property of the digamma function, 



= - 7T cot vr^ - ^ 



(40) 



in order to be able to analytically-continue the OPE to the Minkowski region. Therefore, in 
our specific toy model, duality violations can be computed analytically. Matching onto (|36p . 
one gets 

A(q 2 ) = — 



■ cot 



7T- 



-q 2 + m\ 



(41) 



Recall that indeed the cotangent encodes information about the poles of the digamma func- 
tion, something that the OPE fails to reproduce. Poles are equidistantly located on the real 
axis and are given by the expression 



7T COt 7T 

A similar thing happens with Eq.([3 
wardly and 

A(q 2 ) 



^ oo 



1 



n=l 



(42) 



The analytic continuation can be done straightfor 



4 N r . 



7T 



Jo (qA ) 



(43) 



3 (4tt) 2 

The seemingly negligible exponential terms in the far Euclidean are now poles sitting at 
zeros of the Bessel Jo function. In the Pade approximant, when TV — > oo the residue TZj\r(q 2 ) 
vanishes and we recover the logarithm. In Migdal's approach, in contrast, the poles do not 
merge into the logarithmic cut. 

Interestingly, both the toy model of section 2 and Migdal's model have 5-dimensional 
duals in the context of the gauge/gravity duality, the so-called soft- wall and hard- wall models, 
respectively. It seems interesting to briefly comment on how what we have discussed carries 
over to the fifth dimension. 



5 Connection with holographic QCD models 

In this section we will pursue the relationship between QCD and holographic QCD duals 
suggested in [16]. If Migdal's original work is related to hard- wall holographic models, there 
should be a similar relationship between our toy model and holographic duals with built-in 
Regge behavior. This relationship was recently found in |10| . 
The 5-dimensional model they started from can be written as 



S = - f dtxdze-^y/g-LTt 

J v %i 



(44) 



where the hatted indices jl, v = 0, 1, 2, 3, 4; the metric is parametrised as 

gppdxPdx* = e 2A ^{r 1 ^dxi x dx v + dz 2 ) , i)^ = diag(-l, 1, 1, 1) , (45) 

and $>(z) stands for the dilaton field. For our purposes, namely the spectrum of vector 
mesons, additional bulk fields in the action (|44p are irrelevant. 
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The eigenvalue equation for the vectorial modes can be expressed as a Sturm-Liouville 
differential equation 

where B(z) = $(z) — A(z). After a change of variable, the previous equation can be recast 
in the Schrodinger-like form 

<h + U(z)\ ^ n {z) = m*tl> n {z), U(z) = (?P\ - . (47) 



In the so-called hard-wall models, A = — log z, such that the metric is AdS, and the dilaton is 
taken to be constant, &(z) = 4>, e < z < A, where z = e and z = A are the four-dimensional 
boundaries. Then the potential is given by 

U{z) = ^ ■ (48) 

The Schrodinger equation has Bessel functions as solutions, and the correlator takes the form 
of Eq. (|34jl . yielding m 2 ~ n 2 . The emergence of the scale A is, in the context of holographic 
duals, due to the presence of the infrared boundary brane, which breaks conformal invariance 
and leads to the terms at the outmost right of Eq. (134j) . Duality violations are therefore 
provided by the infrared boundary conditions, and its pattern as a finite-size singularity has 
now a transparent geometrical interpretation: it has to do with the finite size of the fifth 
dimension. 

However, with the geometrical insight of the holographic models we can go further. The 
fact that the eigenfunctions are the Bessel Jo function means that their overlap is non-zero. 
In the language of holography, this leads to non- vanishing decay widths for the resonances. 
Therefore, back to 4 dimensions, the correlated limit of Migdal, which swallows the degree 

of the Legendre polynomials into the argument of the Bessel Jo function, breaks down 
orthogonality and gives an instanton-like size to the resonances. 

All these shortcomings disappear if there is no infrared boundary brane. Then there 
is no instanton-like duality breaking. Moreover, orthogonality of the wave-functions is left 
unbroken and resonances are infinitely narrow, as one expects from a large- N c model of 
resonances. Furthermore, if one modifies the infrared behavior of the model without an 
infrared brane, one can reproduce Regge scaling. This last point was the original motivation 
of the authors of [10J. According to them, the requirement of Regge behavior leads to a 
(asymptotically) unique choice of the dilaton background, namely Q(z) ~ cz 2 . In particular, 
the simplest potential incorporating the right leading order QCD short distances and Regge 
behavior takes the form 

U(z) = c 2 z 2 + JL , (49) 

which corresponds to the three-dimensional harmonic oscillator in quantum mechanics. The 
harmonic term is the dilaton contribution, whereas the centrifugal barrier term is due to the 
AdS metric. The Schrodinger equation is solved through Laguerre polynomials |10j 

^(z)~e-^ 2 z 3 / 2 Li(cz 2 ), (50) 

with the spectrum scaling as 

m 2 n = Acn + Ac, n = 0, 1, ... (51) 
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Furthermore, the decay constants can be computed to give a constant, 



F 2 = \c^F 2 . (52) 
#5 

As noted by |10| . comparison of Eqs. (|5ip and (|52|) with Eq. © shows that the 5-dimensional 
model with linear confinement is dual to the toy model of section 2 with a = 4c and rriy = 4c. 
This immediately means that the dilaton quadratic profile is the holographic dual to the 
duality violating piece 

-q 2 + 4c" 



7T 



A^H^cot 



7T- 



4c 



(53) 



From a phenomenological point of view, this soft-wall model has two free parameters, 55 
and c. 55 can be determined by matching the model to the parton model logarithm, while c 
can be extracted by imposing that rriy = m 2 . Taking m p = 776 MeV, we find 

9 12-7T 2 m 2 9 

g\ = — — , c= — £ ~ 150 MeV 2 . (54) 

It is interesting to notice that the ratio rriy /a = 1 is fixed in the 5-dimensional soft- wall 
model, as opposed to our toy model, where a and my are free parameters. In particular this 
means that the soft-wall model gives a non-zero negative dimension-two condensate. Using 
the OPE expansion of Eq. (|13p one finds 



The model also underestimates the Regge slope, giving a = 0.6 GeV 2 while the phenomenolog- 
ical one hovers around a ~ 1.2 GeV 2 . Finally, the decay constant is predicted to be F = 123 
MeV, while its experimental value lies at F = 210 MeV. Interestingly, both the Regge slope 
and the decay constant adjust to their experimental values with the choice rriy /a = 2. Notice 
that this ratio also cancels the dimension-two condensate. 

It would therefore be really interesting to know if it is possible to modify the soft-wall 
model so that rriy /a = 2. 



6 Conclusions 

Upon reassessing Migdal's work we have shown that, even though the starting point is a 
well-defined Pade approximant to the parton model logarithm, the correlated limits taken 
thereafter not only break the convergence of the approximant but also modify the analytic 
structure of the correlator. We have seen that the correlated limits, which lead to the even- 
tual emergence of an infrared scale A, are nothing but a modeling of quark-hadron duality 
violations. However, this modeling corresponds not to what one expects from a model of 
resonances in large- N c but to finite-size effects, most commonly associated to instantons. Ac- 
tually, the nature of the resonances resembles more that of instanton-like excitations rather 
than quark-antiquark bound states. 

It is worth stressing that in Migdal's approach it is the duality violating piece what brings 
information on the spectrum. This is a general feature: the OPE, being a regular expansion, 
is completely blind to the singularities in the q 2 plane, and therefore unable to reproduce the 
details of the spectrum. Those singularities are encoded in the form of duality breakdown. 
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Therefore, contrary to what was conveyed in Migdal's work, the addition of ultraviolet cor- 
rections to the correlator does not help in determining the spectrum. This is related to the 
well-known fact that different functions can lead to the same asymptotic behavior. As an 
academic exercise, consider our toy model of large- A^ c QCD with built-in Regge behavior. We 
can fine-tune it such that it matches the first k OPE vacuum condensates o Let us imagine 
that all OPE condensates can be accounted for this way. Therefore, the Pade approximant 
(and Migdal's construction upon it) is now unable to distinguish our toy model from real 
QCD. From perturbative QCD and, most generally, the OPE, the hadronic spectrum cannot 
be reconstructed. This same statement was already made in a slightly different context in 
[29]. 

Migdal's original motivation for taking the correlated continuum limit was to make quark- 
hadron duality more accurate. We have shown that it is precisely the other way round: 
quark-hadron duality was enforced from the beginning and preserved by the Pade approx- 
imant. Duality breakdown only occurred after the correlated limits were taken, generating 
the exponentially suppressed terms of Eq. (|35l) . 

Even from a conceptual point of view, in real QCD duality is expected to be restored at 
high energies rather than at low energies. In N c = 3 QCD, and even at large but finite N c , 
resonances in the spectrum develop a width, expected to be proportional to 

IV(n) ~ ^ , (56) 

n being the excitation number. Therefore, increasing widths of resonances makes them 
overlap and eventually the spectrum is smeared out into a continuum: local duality sets in. 
This means that at sufficiently large values of momentum t on Minkowski axis the duality 
picture has to be very accurate^ Actually, based on a model, the authors of [30] found 



Im A(t) ~ exp f - cos (3t (57) 

as a particular realization of Eq. (l37p . a and (3 being constants. However, in the strict large- 
N c limit N c is sent to infinity from the very beginning. Resonances therefore stay infinitely 
narrow no matter the excitation number and local duality never sets in. 

The large- N c and low-q 2 limits taken by Migdal avoid the continuum region of QCD, 
precisely where the assumption of duality is allowed. At low energies resonance poles are 
well-defined only because duality violations are sizeable. In large- N c duality breaks down by 
construction. If we attempt to reconstruct the spectrum of QCD, the assumption of duality 
in Migdal's work is inherently inconsistent. 

The only consistent way to make contact with hadronic resonances with a Pade approxi- 
mant would be to incorporate duality-violating terms in the spectral function before building 
the Pade approximant. Notice, however, that this does not yield any predictions: knowledge 
of the spectrum is already required prior to the Pade construction. Since duality breakdown 
has to be inserted by hand, in model building it is much better to start from an ansatz for the 
spectral function and build the correlator using dispersion relations [21] . if we are working in 
4 dimensions from the start, or playing with the infrared properties of an holographic model, 
if we use 5-dimensional models in AdS backgrounds. 



13 This can be easily achieved, e.g., by relaxing our Eq.© so that the first k decay constants F n are left as 
free parameters. 

14 This is the rationale behind the sum rule program. 
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Gauge/gravity duals provide a more intuitive physical picture of quark-hadron duality 
breakdown, based on geometry. The somewhat misterious infrared scale in Migdal's model 
can now be interpreted as the infrared boundary brane. The presence of this brane breaks 
conformal invariance and generates the terms that we identified as quark-hadron duality vio- 
lations. Both Migdal's model and the hard-wall model start from the parton model logarithm 
and model duality breakdown in exactly the same way. Therefore, the shortcomings men- 
tioned for Migdal's model carry over to the hard-wall model in 5 dimensions. However, in 
gauge/gravity duals there is a way to get rid of those shortcomings by removing the infrared 
brane and changing the infrared properties of the metric and/or the dilaton field at the same 
time. This corresponds to a change in the modeling of the duality breakdown terms. If this is 
done smartly, one can end up with a model with built-in Regge trajectories [10j . A quadratic 
dilaton background then leads to a dual version of our toy model of section 2. Unfortunately, 
the fact that the ratio between m p and the Regge slope in the soft-wall model are fixed leads 
to poor phenomenological predictions. 

The statement that holographic models were leading in an essential way to non-Regge 
behavior [11] was based on the belief that the metric in the ultraviolet was determining the 
spectrum. This is essentially the same philosophy lying behind the approach developed by 
Migdal. In holographic models, as in 4-dimensions, this is not true. Notice that one could 
get rid of the centrifugal barrier in Eq. (|49j) and work in flat space-time with the same dilaton 
background. By solving the corresponding Schrodinger equation the eigenfunctions would 
be proportional to Hermite polynomials, and its scaling would still be Regge-like, m\ ~ n. 
Obviously, the interplay of the metric and the dilaton would give rise to a different background 
and this will affect, e.g., the determination of decay constants and the slope of the spectrum 
(the a parameter in our toy model) but, importantly, will not change its scaling properties. 
Therefore, it is not the ultraviolet of the theory but some modeling of quark-hadron duality 
breakdown what is necessary to model the spectrum. 

This brings up the interesting question of whether there is hope in the AdS-QCD models 
to build the holographic dual of QCD. Unfortunately, very little is known about the form 
of duality violations in real QCDt!£| In order to achieve the goal successfully it seems that a 
better knowledge of non-perturbative QCD is mandatory. 

Acknowledgements 

I am grateful to M. Piai, D. T. Son, M. Unsal, A. Karch, V. Sanz and M. Perez- Victoria for 
useful conversations and to C. P. Herzog, S. Peris and M. J. Strassler for their critical read- 
ing of the manuscript and the suggestions that followed. I also would like to acknowledge 
very fruitful correspondence with S. J. Brodsky, G. F. de Teramond, J. Hirn, G. Lopez- 
Lagomasino, I. Low and E. J. Weniger. My thanks to the Department of Physics of the 
University of Washington at Seattle, where this work has been carried out, for their hos- 
pitality. The author is supported by the Fulbright Program and the Spanish Ministry of 
Education and Science under grant no. FU2005-0791. 

15 See |19l 1211 130] for some recent studies in the large- N c limit. 



13 



Appendix 



As stated in the main text, the vanishing of contour integrals in two-point functions can be 
inferred from the asymptotic behavior of the spectral function, leading in our case to 



1. 



lim -ImnV(t)Qjv(*) ~ f ■ 

t— >oo 7T 



(58) 



However, in the large-iV c limit the above equation is strictly speaking ill-defined. Therefore, 
one has to adopt a slightly different strategy. We start from 

dq 2 



d n 



d(q 



2\n 



n V (q 2 )QN(q 2 ) 



ill 
2tt7 

n\ 
n\ 



(q 2 + n 2 )^ 1 
dq 2 

, g2|=S0 (^+ M 2)n + 1 



n V (q 2 )Q N {q 2 ) 

n v (q 2 )Q N (q 2 } 



+ 



no 



dt 



o 



(t + n 



2\n+l 



IT 



lmU v (t)Q N (t) 



(59) 



where the limit sq — ► oo is to be understood henceforth. To get to (|26|) in the main text we 
have to show that the integral in the second line above identically vanishes. Even though we 
will carry all explicit calculations with the model of section 1, we will find that this statement 
is a general result of large- N c QCD itself. 
Our model reads 



Kv(q 2 ) = — 



-q 2 + m\ 



-Im n y (f) = F 2 \^ S(t — my — an) 

7T — i 



(60) 



n=0 



Recall that the digamma function VKO i s analytic everywhere in the complex £ plane except 
for simple poles at negative integer values of the argument, i.e., £ = 0, —1, —2, .... Moreover, 
it satisfies 

Resfo&(£);£] = -1, £ = 0,-1,-2,... (61) 

The second line in (|59p can then be straightforwardly evaluated by reiterative application of 
Cauchy's residue theorem. The integrand has a pole of order n+1 in the Euclidean axis sitting 
at q 2 = —n 2 along with an infinite number of equidistant (simple) poles in the Minkowski 

cLXIS , 3/t 

q 2 = m 2 v + ak, = 0, 1, ... (62) 

Application of the residue theorem then yields 
dq 2 



ill 



2mT lq2l=S0 {q 2 + a 2 Y+^ 



n V {q 2 )Q N {q 2 ) 



+ lim 

q 2 — *m^-\-ak 



lim 



d" 



M 2 d{q 



2\n 



Kv{q 2 )Q N (q 2 ) 



+ 



A:=0 



2 - (my + ak) 
(q 2 + i^+I 



Uy( q 2 )Q N (q 2 



. (63) 



Obviously, the first term is the left-hand side we started from in (|59p . Therefore, the second 
line above has to cancel the integral over the spectral function in (|59p . Indeed, 



*0 



dt 



(t + » 2 ) 



2\n+l 



-ImU v (t)Q N (t) 

TT 



QN(m v +ak) 



k=0 



(my +ak + /i 2 ) n+1 



(64) 
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whereas 



lim 

g 2 — *my+ak 



E 

k=0 



(m v + ak) 



(q 2 + [i 2 



<-^Il v ( q 2 )Q N { q 2 



lim F 2 V 

2 i „L. 



k=0 



-q z +m v 



+ k)i/> 



-q 2 + m 2 v 



Qn(q 2 



2\n+l 



(65) 



which equals (|64p but for a change in sign] 16 ! 

Now we only need to show that the first line in (|63[) cancels (|65p . or equivalently, as we 
have just checked, equals (|64p . In other words, the contribution from the infinite poles on the 
Minkowski axis counterplays that of the pole at the Euclidean axis. Applying the binomial 
formula to the first line of (1631) . one finds 



d(q 



2\n 



F" 



-q 2 + my 



F 2 



E 

~ 3=1 

oo n 

f2 EE 



l\3 



-q 2 + m v 



n 



j)\ (m 2 y + ak + fi 2 )3+i ' 



(66) 



k=0 j=l 

where in the last line we have used the definition of the polygamma function, namely 

k=0 



' ^2, 1 



(67) 



The last line in (|66[) can be recast in the more useful way 



E 



(my + ak + [i 2 



k=o ^-°y 

Therefore, in order to match 



\n+l E 



nl 



\ ( n ~ j) 



.Q%->\-V 2 )(m 2 v + ak + n 2 ) n - 1 



it follows that 



Qjv( m 



v 



n 1 



-Q 



(n-j) i 



-H 2 )(m v + ak + fi 



2\n-j 



(68) 



(69) 



but the sum over j is just the Taylor expansion of Q,N{q 2 ) around q 2 = — /i 2 , evaluated at 
q 2 = my + an, and truncated after n + 1 terms. However, being Q.N{q 2 ) a polynomial of 
degree N, the previous relation is exact if n > N, which is one of our starting points. This 
completes the proof. 

Recall that the proof, as anticipated, relies only upon the analytic properties of two- 
point Green's functions in large--/V c , namely the existence of no cuts and just simple poles. 
Therefore, the result is model- independent and valid for large- iV c QCD. 



1 This becomes obvious if one deforms the path 7 such it does not include the pole in the Euclidean. Then 
Cauchy's integral theorem asserts that the integral over the circle equals that over the physical axis. 
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